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Abstract 

We study the cosmology of a galileon scalar-tensor theory, obtained by covariantizing the 
decoupling lagrangian of the Dvali-Gabadadze-Poratti (DGP) model. Despite being local in 
3+1 dimensions, the resulting cosmological evolution is remarkably similar to that of the full 
4+1-dimensional DGP framework, both for the expansion history and the evolution of den- 
sity perturbations. As in the DGP model, the covariant galileon theory yields two branches 
of solutions, depending on the sign of the galileon velocity. Perturbations are stable on 
one branch and ghost-like on the other. An interesting effect uncovered in our analysis is 
a cosmological version of the Vainshtein screening mechanism: at early times, the galileon 
dynamics are dominated by self-interaction terms, resulting in its energy density being sup- 
pressed compared to matter or radiation; once the matter density has redshifted sufficiently, 
the galileon becomes an important component of the energy density and contributes to dark 
energy. We estimate conservatively that the resulting expansion history is consistent with the 
observed late-time cosmology, provided that the scale of modification satisfies r c > 15 Gpc. 
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1 Introduction 



Scalar-tensor theories of gravity have experienced a resurgence of sorts, over the last twenty years. 
This is due in part to string theory, where the plethora of compactification moduli generically 
appear in the 4D effective theory with kinetic mixing with the graviton. Moreover, the discovery 
of accelerated expansion makes the possibility that General Relativity is modified on the largest 
scales plausible. If this is the case, then the new gravitational degrees of freedom relevant on 
cosmological scales are likely to include a scalar cousin for the graviton. 

The best-known example of a scalar-tensor theory is due to Brans and Dicke (BD) [1], 

Sbb = ^J d'xV^ (*R - (d$) 2 ) + J d 4 xV^~gC mattcr [g] , (1.1) 

where the matter Lagrangian is independent of $. Unfortunately, the BD parameter is so tightly 
constrained by solar system and pulsar observations, u;bd ^ 4 x 10 4 [2], that the cosmological 
effects of the BD scalar are rendered uninterestingly small. 

A tantalizing alternative is that the apparent decoupling of the scalar field is a local effect, 
owing to the large matter density of the solar system or pulsar environment. In other words, the 
BD parameter is effectively a growing function of the density. While decoupled locally, the scalar 
field can have interesting cosmological effects in the much sparser cosmic environment. There are 
only two robust mechanisms that realize this idea. One is the chameleon mechanism [3, 4, 5]: by 
adding a suitable potential V(<&), the scalar field acquires mass which depends on the density. The 
mass is large in regions of high density, thereby suppressing any long-range interactions. Theories 
of f(R) gravity [6] rely on the chameleon effect to ensure consistency with solar system tests [7]. 

A second mechanism is the Vainshtein screening effect of the longitudinal graviton or brane- 
bending mode, usually denoted by ir, in the DGP model [8]. As we review in Sec. 2, this effect 
is most easily understood in a certain decoupling limit of the theory [9, 10]: M P \,M 5 — > oo, 
keeping the strong coupling scale (Mpir" 2 ) 1 / 3 fixed. The resulting theory is local on the brane, 
and describes a self-interacting scalar field coupled to weak-field gravity in 3+1 dimensions: 

M 2 r 2 1 

C = --^h^iSh)^ + Ml^(Eh), u - -^-(dn) 2 Dir + -hTT^ , (1.2) 

where £^h a p = + . . . is the linearized Einstein tensor. As a vestige of 5D Lorentz 

transformations, the 7r action is invariant under the Galilean shift symmetry, d^n — > <9 M 7r+c M . Thus 
7r has been dubbed a galileon field [11]. In regions of high density, p ^> Mp Y r~ 2 , non-linearities in 7r 
dominate and result in its decoupling. This is qualitatively similar to the chameleon mechanism, 
except that the galileon relies on derivative interactions as opposed to a scalar potential. 
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In this paper, we study the cosmology of the galileon, by promoting (1.2) to a fully covariant, 
non-linear theory of gravity coupled to a galileon field: 



S = J A'x^g (|i e - 2 *fl - {^(dn) 2 nn + C^fo]) . (1.3) 

As in the DGP model, where the Galilean shift symmetry is only exact in the strict decoupling 
limit Mpi — > oo, the shift symmetry is now broken by M P1 -suppressed operators in (1.3). Of course 
the above non-linear completion is by no means unique — many other Lagrangians, for instance 
including a (97r) 4 /Mpj term, will reduce to (1.2) in the weak-field limit. Since the cosmological 
predictions should be fairly robust under such corrections, however, we take (1.3) as a fiducial 
galileon theory and study its implications for cosmology. 

Despite being a local theory in 3+1 dimensions, the cosmology derived from (1.3) comes 
remarkably close to reproducing that of the 4+1-dimensional DGP model, at least for Hr c > 1. 
As we will see, the agreement holds for both the expansion history and the evolution of density 
perturbations. For a sneak preview, see Figs. 6 and 7, respectively. More generally, our galileon 
cosmology reproduces many qualitative features of DGP: 

1. The Friedmann equation allows for two branches of solutions, depending on the sign of n. 
In analogy, the modified Friedmann equation in DGP [12] also has two branches: 

2. One branch of solutions has stable perturbations, whereas the other is plagued with ghost- 
like instabilities. This again agrees with DGP [10], where the "— " and "+" branches in (1.4) 
are stable and unstable, respectively. 

3. The effective equation of state for the galileon satisfies w n < —1 on the stable branch, and 
w n > — 1 on the unstable branch. This agrees with the effective equation of state inferred 
from the H/r c correction in (1.4) [13]. 

4. Moreover, the two branches are classically disconnected, unless R < 0. This is closely related 
to the condition p — 3P < — YlM^/r^ necessary to transition from the stable to the unstable 
branch of solutions in the decoupling theory [10]. 

An important difference with DGP, however, is that our covariant galileon theory does not allow 
for self-accelerated cosmology — the self-accelerated solution is spoiled by l/M P \ terms in (1.3). 

One of our key results is a cosmological analogue of the Vainshtein screening mechanism. At 
early times, Hr c ^> 1, the dynamics of n are dominated by the cubic interaction term, resulting 



in the galileon energy density being suppressed by 0(1 /Hr c ) compared to the matter or radiation 
fluid. When the matter density has dropped sufficiently, so that Hr c ~ 1, the galileon becomes 
an important component of the total energy density and contributes to dark energy. 

We also study the effects of the galileon on the growth of inhomogeneities. By virtue of its non- 
minimal coupling to gravity, the galileon enhances the gravitational attraction between particles, 
which translates into more efficient growth of density perturbations. The screening mechanism 
is also at play in the evolution of perturbations: the galileon enhancement is suppressed for 
Hr c ^> 1, but becomes important once Hr c ~ 1. A similar time-like Vainshtein effect was also 
observed in [14]. 

While a full likelihood comparison to data is left for future study, we discuss various constraints 
on the galileon cosmology, such as from estimates of the matter density at different redshifts, the 
luminosity distance relation, and the angular-diameter distance to the last scattering surface. The 
resulting bound on r c is 

r c >15Gpc, (1.5) 

which constrains the scale of the modification to be at least a few times the Hubble radius today. 
However, since all other cosmo logical parameters are kept fixed in our considerations, (1.5) is likely 
a conservative estimate. 

The paper is organized as follows. In Sec. 2 we review how the weak-field action (1.2) arises from 
the decoupling limit of the DGP model, and describe the origin of the self-screening mechanism 
near spherical sources. In Sec. 3 we discuss the non-linear extension (1.3) and derive the covariant 
equations of motion. We present in Sec. 4 the cosmology of the galileon model. In particular, we 
derive an approximate analytic solution, which displays the screening mechanism, and study its 
stability. Our analytic considerations are borne out by the numerical solutions presented in Sec. 5. 
Turning to inhomogeneities, we study in Sec. 6 the effects of the galileon on the growth of density 
perturbations. In Sec. 7 we discuss various observational constraints on galileon cosmology and 
derive the bound on r c given in (1.5). We conclude in Sec. 8 with a brief summary and discuss 
future research avenues. 

While preparing this manuscript we became aware that Deffayet et al. were independently 
studying a model with some similarities to ours [15, 16]. This paper is based on the M.Sc. thesis 
of N.C. at the University of Waterloo [17]. 

2 Decoupling Limit of DGP 

In the DGP model, our visible universe is confined to a 3-brane in a 4+1-dimensional bulk. Despite 
the fact that the extra dimension is infinite in extent, 4D gravity is nevertheless recovered over 
some range of scales on the brane because of a 3+1-dimensional Einstein-Hilbert term, intrinsic 
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to the brane: 



r m 3 f [ M 2 \ 

S DGP = / tfxy/^-f-Rs + / d*xy/=& ~^Ra + C matter [g] . (2.1) 

J bulk 1 Aranc V 1 J 



The bulk and brane Planck masses define a cross-over scale, 

which separates the AD and 5D regimes. At distances r C r c on the brane, the gravitational force 
law scales as 1/r 2 , whereas for r ^> r c it scales as 1/r 3 . 

From the point of view of a brane observer, the 5 helicity-2 states of the massless 5D graviton 
combine to form a massive spin-2 representation in AD. More precisely, the AD graviton is a 
resonance — a continuum of massive states — whose spectral width is peaked at the scale r~ l . 
As in massive gravity [18], the helicity-0 or longitudinal mode, denoted by n, becomes strongly 
coupled at a much lower scale than M 5 , given by [9, 19] 

A strong = (M P1 r c ~ 2 ) 1/3 . (2.3) 

Forr c ~ Hq 1 = 10 28 cm, for instance, this gives A st J ong ~ 1000 km. The strong coupling behavior is 
essential to the phenomeno logical viability of the model through the Vainshtein screening effect [18, 
19]. As we review below, non-linear interactions in n are important near an astrophysical source 
and result in the decoupling of rr from the source. The characteristic scale below which ir is 
strongly coupled, denoted by r*, is given by 

r* = (r c 2 r Sch ) 1/3 , (2.4) 

where r^ch is the Schwarzschild radius of the source. And since r c is cosmologically large (of order 
of the Hubble radius today), r\ is parametrically larger than r Sch . 

In analogy with massive gravity [20], it is instructive to zoom in on the non-linearities in 7r by 
considering the decoupling limit [9, 10]: M P i, M 5 — > oo keeping the strong coupling scale A str0 ng 
fixed. Equivalently, around a spherical source this corresponds to sending rs c h — ► keeping 
fixed. In other words, in this limit non-linearities in the helicity-2 (Einsteinian) modes drop out, 
while interactions of the helicity-0 state survive. The resulting effective theory is local on the 
brane and describes (weak-field) gravity plus a scalar field n in 3+1 dimensions: 

M 2 ~ r 2 1 ~ 1 

^Einstein = ~ ~^ ^ (£ h) , u - 3(9tt) 2 - ^(^DtT + -h^T, u + — TtT , (2.5) 

where hap — + ... is the linearized Einstein tensor. This lagrangian is, up to a total 

derivative term, invariant under the Galilean shift symmetry, 

<9 M 7r -> <9 m tt + c M , (2.6) 

which is a vestige of the full 5D Lorentz transformations. Thus n has been dubbed a galileon 
field [11]. 
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2.1 Self-screening effect 



The approximate recovery of general relativity in the vicinity of astrophysical sources, through 
the Vainshtein effect, can be understood at the level of (2.5) [10]. The equation of motion for the 
galileon, 

9" (6M P1 ^vr + 2r c 2 ^7rD7r - r^dvr) 2 ) = —T , (2.7) 

is remarkable in many respects. Even though the interaction term in (2.5) contains four derivatives, 
the equation of motion is nevertheless second-order — all higher-derivative terms cancel out when 
performing the variation. Moreover, (2.7) takes the simple form d^j^ = —T/2Mp\ for some 
7r-current jj*, thereby allowing for a generalized "Gauss' law": spherically-symmetric exterior 
solutions for n only depend on the mass enclosed. 

Let us indeed study the spherically-symmetric galileon profile due to a point mass: T = 
—M5 3 (r). In this case, (2.7) can be integrated to give 

6tf,^(r)+*^ = ^. (2.8) 
Using (2.4), the solution to this algebraic equation for ir' is given by 

7r'(r) 



Mi 



pi 




(2.9) 



Note that we have chosen the branch of the solution such that it' — > as r — > oo. The other 
branch, corresponding to n' diverging at infinity, belongs to the same branch of solutions as 
the self-accelerated DGP cosmology, and is therefore unstable [10]. This is a general property: 
solutions to (2.7) always come in a pair, with one member continuously connected to the trivial 
solution, with stable perturbations, and the other connected to the self-accelerated cosmological 
solution, with unstable perturbations. It is impossible to classically move from one branch of 
solutions to the other without violating some energy condition [10] . In this work, we focus almost 
exclusively on the stable branch of solutions. 

At short distances, r <C r*, the galileon-mediated force is clearly suppressed compared to the 
gravitational force: 



7-i I 3/2 2 / \ 3/2 

F n V7T ry r I r N 



^grav Afpi|V$| r c rl/2 



<1. (2.10) 



Thus, as advocated, the strong interactions of 7r lead to its decoupling near a source, and the theory 
reduces to Newtonian gravity. This approximate recovery of standard gravity near a source has 
been established in approximate solutions of the full DGP model [19, 21, 22]. The above ir- 
mediated force, albeit small in the solar system, is nevertheless constrained by lunar laser ranging 



observations [23, 24, 25]: r c > 120 Mpc. A comparable bound on r c has also been obtained by 
studying the effect on planetary orbits [26] . 

At large distances, r ^> r*, on the other hand, the non-linear terms in n are negligible, and 
the resulting correction to Newtonian gravity is of order unity: 



F 1 

T^ = \- ( 2 - n ) 

r grav <J 

The galileon-mediated force therefore leads to an enhancement of the gravitational attraction by 
a factor of 4/3. In this far-field regime, the theory reduces to a scalar-tensor theory, with the 
galileon acting as a Brans-Dicke scalar. 

2.2 Jordan frame description 

Our action (2.5) is cast in Einstein frame, where the kinetic terms are diagonal, but n couples 
directly to matter. We find it more convenient to instead work in Jordan frame, by performing 
the shift 

2,71 

h^ u = h^ v + —— Tj^i/ . (2.12) 

This removes the nT coupling, at the price of introducing kinetic mixing between h and n: 

M 2 r 2 1 

Jordan = --^h^(£h)„ u + M^^Sh)^ - -f-(dn) 2 Dir + -bTT^ . (2.13) 

This forms makes the Brans-Dicke nature of the theory manifest, in the limit where the 7r- 
interactions can be neglected, with the Brans-Dicke parameter identified as c^bd = 0. 



3 Non-linear Completion 

Nearly all of the interesting phenomenological features of the DGP model are attributable to the 
helicity-0 mode tt and can be understood at the level of the decoupling theory. The Vainshtein 
effect, reviewed above, is one example. The existence of a self-accelerated solution is another 
example: the equation of motion (2.7) in vacuum (T = 0) has a solution where 7r ~ M-pix^x^ /r 2 , 
in agreement with the weak- field limit of de Sitter space [10]. 

This motivates us to propose a AD theory of modified gravity, by promoting (2.13) into a fully 
covariant, non-linear theory of gravity coupled to a galileon field. By construction, this non-linear 
theory will reduce to (2.13) in the limit of weak gravitational fields. Therefore its predictions will 
agree with those of the full DGP model to leading order in 1/M P1 . 

Looking at (2.13), a natural non-linear completion suggests itself: 

S = J d'x^-g (Ml e -WM PljR _ ^-(dnfOn + £ mattcr [^ , (3.1) 
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where □ is now understood as the covariant Laplancian: □ = V M V M . This clearly reduces to (2.13) 
in the weak-field limit. The Galilean shift symmetry (2.6) is softly broken in the action (3.1) 
through Mpi-suppressed operators. This also true of the full DGP model, where the Galilean 
symmetry arises only in the strict decoupling limit as a remnant of the full 5D Lorentz group. 

Of course the above non-linear extension is by no means unique. For instance, we could 
consider more general functions of n/Mp\ multiplying the Ricci scalar, or include corrections of 
the form (d7r) 4 /Mpj, all of which would drop out in the limit Mp\ — > oo. Be that as it may, we 
take (3.1) as a fiducial covariant theory and explore its cosmological predictions. A study of more 
general lagrangians is left for the future. 

Remarkably, as we will see in the next Section, the AD cosmology arising from (3.1) reproduces 
many features of the full-fledged DGP model. Our Friedmann equation has two branches of 
solutions, depending on the sign of the velocity of n. The two branches are distinguished by having 
stable or unstable (ghost-like) perturbations. Moreover, we uncover a cosmological analogue of the 
Vainshtein effect: at early times, when the density of the universe is high, non-linear interactions 
in it are important, resulting in the galileon energy density being subdominant compared to the 
matter or radiation fluid. 

The covariant equation of motion for the galileon is readily obtained from (3.1): 

M 2 

(□tt) 2 - (V m V,tt) 2 - it^V^V.TT = -^Re' 2 ^ . (3.2) 

Similarly, the Einstein equations are given by 

e- 2 ^>M^(? F = T^ + M 2 1 (V^-^D)e-^ /Mpi 

2 

+ -^(2V M 7rV,7rD7r + ^V Q 7rV a (97r) 2 -2V (M 7rV,)(<97r) 2 ) . (3.3) 
Mpi 

Since the matter action is independent of ir, the matter stress-energy tensor satisfies the usual 
conservation law: V M T U „ = 0. 



4 Cosmology 

In this Section we specialize the above equations to the cosmological context, by assuming homo- 
geneity and isotropy. For simplicity, we focus on the case of a spatially-flat universe. Under these 
assumptions, (3.2) reduces to 

^ (Hit 2 ) + 3H (Hit 2 ) = ^-Re- 2 ^ , (4.1) 
at mf, 

where dots represent derivatives with respect to proper time t. Remarkably, the left-hand side is 
reminiscent of the equation for a canonically-normalized scalar field, + 3H<p, with <fi identified 
as Hit 2 . 
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Because the effective field momentum is proportional to 7r 2 , however, the dynamics are quite 
different from those of a standard scalar field. In particular, if R > 0, as is the case in a universe 
dominated by matter, radiation or vacuum energy, then solutions with 7T > and tt < are 
classically disconnected. Indeed, expanding (4.1), 

M 2 

2Hnn + ... = ^R e - 2 ^ , (4.2) 
6r 2 

we see that -k is driven away from zero if R > 0. A necessary condition to transition from one 
branch to the other is therefore R < 0. This is closely related to the condition p — 3P < — YlM^/r 2 
necessary to transition from the stable to the unstable branch of solutions in the decoupling 
theory [10]. We will see that the above conclusions are borne out by numerical analysis — as long 
as the universe is dominated by matter, radiation, or vacuum energy, n never changes sign. 

We next turn our attention to the Friedmann and Raychaudhuri equations for the scale factor. 
For the matter, we assume as usual that the stress-energy tensor is described by a perfect fluid, 

= {p + P)u^u u + Pg^ , (4.3) 

with p and P denoting the energy density and pressure of the fluid. The Friedmann equation is 
given by the (0,0) component of (3.3): 

SM^tfV 2 ^ 1 = p + QHttM pi (^e- 2 ^ Mpi - r -^Pj . (4.4) 

Similarly, the Raychaudhuri equation follows as usual from a combination of the (i, i) and (0, 0) 
components: 

M 2 ie - 27 r/M P1 « = 
(X 

+ 

From the form of (4.4) and (4.5), we can read off an effective energy density and pressure for 
the 7r field, which we denote by p w and P n . Of course, since n is non-minimally coupled to gravity, 
p n is not conserved, hence it should only be understood as an effective energy density, informing 
us about the effects of the galileon on the cosmological evolution. In any event, we find: 

P7T = 6HirM Pl (V 2 ^ - ^ ; 

P w = 2tt 2 f 7 ^- + 2e- 2 ^ M A - 2M P1 e- 2 ^ Mpi (tt + 3#tt) 
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M 



I>1 



(4.5) 



(4.6) 



Note that p n changes sign depending on the choice of branch, i.e., whether tt is positive or negative. 
That the galileon effective energy density can be negative should not come as a surprise, since it 
is well-known that non-minimal couplings can induce violations of various energy conditions in 
Jordan frame [27, 28, 29]. This can also be seen at the level of the effective equation of state, 



which, a priori, allows for < — 1. In fact, we will see in Sec. 4.2 that w w < —1 holds at early 
times on the stable branch, when the universe is radiation- or matter-dominated. We will also see 
in Sec. 5.1 that w n anyway contributes positively to the total effective equation of state, w to t (see 
Fig. 4), since p w is negative. 

Many of these features also arise in the full DGP model. Indeed, if we think of the H/r c 
modification in the DGP Friedmann equation [12], 



as an effective contribution to the matter content, then clearly its energy density can have either 
sign, depending on the choice of branch. By the same token, the effective equation of state 
corresponding to the modification can be < —1. For instance, on the normal branch, the effective 
energy density is negative and w e s < — 1 [13]. While such phantom behavior may at first seem 
surprising, our analysis now makes it clear that it is nothing but a natural consequence of the 
scalar-tensor nature of gravity on the brane. 

4.1 Self-accelerated solution? 

In analogy with the full-fledged DGP model, we are tempted to look for a self-accelerated cos- 
mology — a de Sitter solution in the absence of any cosmological constant or matter other than 
7r itself. In the regime \tt\ <C Mp\ where we expect agreement with DGP, the equations of motion 
at first sight do seem to allow for approximate self-acceleration. Setting -k = 7r and H = H to 
be constant in (4.1) and (4.4), one finds an approximate solution for \ir\ <C M Ph given by 



where the "?«"'s indicate the assumption \n\ <C Mp\. Up to a trivial redefinition of r c , this agrees 
with the self-accelerated DGP solution that follows from the "+" branch of (4.8): Hq GP = l/r c . 
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(4.7) 




(4.9) 



Unfortunately, the approximation | tt | Mpi breaks down within a time t ~ r c , which is of the 
order of a Hubble time. In other words, H evolves significantly over a Hubble time, and hence 
cannot be approximated as constant. Self-accelerated cosmology is spoiled in our galileon theory 
by n/Mpi corrections. 

4.2 Early-time solution and cosmological screening 

In this Section we derive approximate analytic solutions for when the universe is dominated by 
other components than n, such as matter, radiation or dark energy. We will see that the dynamics 
of the galileon exhibit a time-like analogue of the Vainshtein effect: at early times, t <C r c , non- 
linearities in n are important, resulting in the galileon energy density being negligible. Once 
t ~ r c , however, n exits the strongly-coupled regime, and becomes a significant contribution to 
the total energy density. 

To simplify the analysis, suppose that the universe is dominated by a single matter component 
with constant equation of state w. Moreover, we assume that the variation in n throughout this 
phase is small in Planck units: |A7r| <C Mp\. Since we can always set ir(t = 0) = by trivial 
rescaling of Mp\, it follows that e 2?r / Mpi ~ 1. The consistency of these approximations will be 
checked a posteriori. 

With these assumptions, the Friedmann equation (4.4) reduces to its standard form, 3if 2 Mp : 1 « 
p, with the usual solution 

a(t)^t^^. (4.10) 

Substituting this into (4.1), 

we see that the galileon equation of motion allows for a solution with n = constant: 

7T 1 1 — 3W 

Mpi = 7 c p(l-w) 

Therefore a real solution exists if either w < 1/3 or w > 1, but, as we will see shortly, only for 
w < 1/3 is the solution a dynamical attractor. Moreover, we will see in Sec. 4.2.3 that among the 
two branches of solutions in (4.12), only n < has stable fluctuations. The other branch, with 
7T > 0, is plagued with ghost-like instabilities. 
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(4.11) 



(4.12) 



4.2.1 Cosmological screening 



The above constant-7r solution is only valid provided that the galileon energy density is a negligible 
contribution to the total energy. With our approximation e 2w ^ Mpi ~ 1, the first of (4.6) reduces to 



Pn 


2 tV 


(- 


r 2 c TX 2 


3H 2 M 2 X 


H M P1 








/- 


- 3w 




± Hr c \ 


'3(1 


— w) 



(4.13) 



This elucidates the time-like screening effect advocated earlier: at early times, when Hr c ^> 1 (or, 
equivalently, t < r c ), non-linear galileon interactions are important, and as a result its gravita- 
tional backreaction is negligible. In particular, this ensures that nucleosynthesis and recombination 
proceed as in standard cosmology, with negligible corrections coming from n. We note in passing 
that the parametric dependence of p^/SH^M^ is consistent with the modified Friedmann equation 
in the full DGP model: looking back at (4.8), the relative contribution of the H/r c modification 
term to the total expansion rate is indeed suppressed by 1/Hr c . 

The self-screening effect breaks down after a time of order r c , at which point p n becomes a 
significant contribution to the expansion rate. As we will see in Sec. 7, constraints on the late-time 
expansion history will enforce a lower bound on r c . Incidentally, t ~ r c also signals the moment 
when the approximation \ir\ <C M P \ breaks down, since 



7T 



Ml 



PI 



l ~^ W * (4.14) 



3(1 — w) r 



c 



To summarize, the constant-^ solution in (4.12) is valid at early times, t r c , when the galileon 
is strongly coupled. In this regime, the galileon energy density can be consistently neglected, and 
the galileon excursion in field space is small in Planck units. 

An illustrative way to understand this cosmological screening is to consider the non-linear 
galileon interactions as an effective BD parameter. By comparing the BD and galileon actions, 
given by (1.1) and (3.1) respectively, and making the identification $ = e~ 27T ^ Mp \ we can define 
an effective dynamical u^, 



1 /3(l-3w), r 

where in the second line we have made the approximation e 27r / Mpi ~ 1, and substituted the 
constant-7r solution (4.12). Following the discussion of the previous paragraph, oj^q is large at 
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early times, when Hr c ^> 1, and standard cosmology is recovered. When Hr c ~ 1, however, we 
enter the scalar-tensor gravity regime with ~ 1. 

For completeness, we also derive the galileon equation of state during the strong coupling 
phase. Substituting (4.12) into (4.7), we obtain 

„. = -?(! -„)+£, (A.). (4.16) 

Interestingly, in this regime w n is completely fixed by the background equation of state. In 
particular, this says that w n < — 1 for w < 1/3. As discussed earlier, however, this does not 
signal the presence of ghost instabilities, but instead is a natural consequence of the non-minimal 
coupling to gravity. 

4.2.2 Dynamical Attractor 

We next prove that our constant-^ solution is a dynamical attractor for physically-relevant values 
of w. Perturb the galileon as 

7T(*)=7f(t) +¥>(<), (4.17) 

with 7r given by (4.12). Since the backreaction of p w on the expansion has already been shown 
negligible in this strongly coupled regime, H{t) and R{t) are oblivious to ip and hence can be left 
unperturbed. Expanding (4.1) to linear order in </?, while remembering that e 2w ^ Mpi ~ 1 in our 
approximation, we get 

{Hip) + 3H {Hip) « . (4. 18) 

It follows that Hip ~ I /a 3 , or 

1 _l = w 

(p ~ — — ~ t n-» . (4.19) 
Ha A 

Thus perturbations redshift away compared to it = constant for — 1 < w < 1, and hence the 
solution is an attractor in this range. On the other hand, from (4.12) we know that the constant-7r 
solution only exists for either tu<l/3oriu>l. The above analysis has now established that the 
solution is stable for w < 1/3 and unstable for w > 1. 

4.2.3 General stability analysis 

We can easily extend the above stability analysis to general perturbations: 

n(x,t)=n(t) + <p(x,t). (4.20) 

We work at the level of the action, since here it is transparent whether perturbations have a right- 
sign kinetic term or are ghost-like — a diagnosis that is trickier to make with linearized equations 
of motion. 
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Since it is strongly coupled at early times, its perturbations are dominated by the cubic inter- 
action term in (3.1): 

SeS = ~itJ d4a V=0(frO 2 a* • ( 4 - 21 ) 

Once again, as in the analysis of Sec. 4.2.2, we neglect perturbations in the metric, since the 
cosmological background is driven by some other source of stress-energy. Expanding (4.21) to 
quadratic order in <p, then after some integration by parts we obtain 

SS = J d A x^ (Dng^ - V m V,tt) d»ipd»<p . (4.22) 

Note that this holds for arbitrary background n(x, t). To make contact with the decoupling results 
of [10], it is useful to define K^ v = — r^V^V^^r/Mpi, which measures the extrinsic curvature of the 
brane in the decoupling limit of the full DGP model. In terms of K^, (4.22) takes the form 

= -2 J d 4 x^g~ (K^ - g^K) WipVip . (4.23) 

In the weak-field limit, where we expect consistency with the decoupling limit, this indeed agrees 
with Eq. (22) of [10] in the strong coupling regime: K^ u ^> 1. Remarkably, the form of the moduli 
space metric, K^ u — g^ u K, is preserved in our AD covariant theory. 
Specializing to the constant-^ solution, (4.22) reduces to 

S$ = ^J d^x^-gHn {-^ + jj(Vv?) 2 ) . (4.24) 

Thus the kinetic term of fluctuations is proportional to 7f — the kinetic term is positive if tt < 
0, corresponding to stable perturbations, and is negative if tt > 0, corresponding to ghost-like 
perturbations. Hence, just like in the full DGP model, one branch of solutions is stable, whereas 
the other is unstable. Moreover, looking back at (4.13), the stable branch has p n < 0, whereas 
the unstable branch has p n > 0. Again this is consistent with DGP — interpreting the H/r c 
modification in the DGP Friedmann equation (4.8) as an effective energy density, then this energy 
density is negative on the stable (minus-sign) branch and positive on the unstable (positive-sign) 
branch. 

5 Numerical Analysis 

We can solve the 7r equation of motion (4.1) and the Raychaudhuri equation (4.5) numerically to 
obtain exact cosmological solutions. We focus exclusively on the stable branch, by setting initial 
conditions with -k < 0. Since p n < on this branch, we must include a dark energy component to 
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1+z 



Figure 1: Result of numerically integrating the cosmological evolution equations with r c = 10 Gpc. 
The solid curves denote the matter, radiation, cosmological constant and galileon energy density. 
The dotted line is the galileon energy density as predicted from the analytic solution. 




1+z 

Figure 2: Same as Fig. 1, this time for r c = 20 Gpc. 
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obtain late-time acceleration, which we take to be a cosmological constant A for simplicity. In 
other words, in addition to -n we include radiation, matter and a cosmological term. And since 
contributes to the effective dark energy component, for each value for r c we adjust A so that the 
fractional contribution in matter today is kept fixed to the fiducial value = 0.26. 



5.1 Results 

Figures 1 and 2 show the evolution of the energy density for the various components as a function 
of redshift, for r c = 10 and 20 Gpc, respectively. (Since these are log-log plots and p w < 0, for 
the galileon component we instead plot \p n \.) These figures confirm the analytic results derived in 
Sec. 4.2. The dotted and dashed lines in Figs. 1 and 2, respectively, trace the galileon energy den- 
sity as predicted by the constant-^ solution of Sec. 4.2. Note that, to plot the analytic prediction, 
we have substituted in (4.12) the total equation of state w tot defined by 

J^ = -|(l + «**)■ (5-1) 

From the figures, we see that for most of the cosmological evolution, the actual galileon energy 
density agrees well with the constant-^ prediction. Moreover, the constant--^ solution is manifestly 
an attractor: as seen in Fig. 2, for instance, the galileon energy density quickly converges to the 
analytic prediction. The exact solution starts to deviate from the analytic prediction around 
the present time, however, when the Hubble radius if -1 becomes comparable to r c — this is 
consistent with the discussion of Sec. 4.2.1. Note that the galileon energy density is a more 
significant contribution to the dark energy today for r c = 10 Gpc than for 20 Gpc. 

Figure 3 shows the fractional contributions to the total energy density, defined for the various 
components as 

2tt/M p1 

= (5 ' 2 » 

Note that the galileon contribution is negative, since tt < 0. This figure shows once again that the 
galileon is subdominant at early times, until Hr c ~ 1 when p n becomes a significant contribution 
to the expansion rate. 

Plotted in Fig. 4 is the effective equation of state for the whole evolution, w tot , defined in (5.1). 
Here we plot the ACDM behavior (solid curve) compared to our galileon model with r c = 10 Gpc 
(dotted curve), 15 Gpc (dash-dotted) and 20 Gpc (dashed). The early-time behavior is as expected 
from standard cosmology — the equation of state goes through successive stages of w to t ~ 1/3 
and Wtot ~ 0, corresponding to radiation- and matter-dominated eras, respectively. At late times, 
the energy density is dominated by an effective dark energy component composed of A and p n , 
with negative equation of state. Compared to ACDM, the galileon contribution pushes the dark 
energy equation of state to values larger than —1. And the smaller r c is, the further iy tot is from 

15 



-0.5 
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Figure 3: The fractional contributions to the total energy density from matter (fi m ), radiation 
(f2 r ), cosmo logical constant (Qa, dashed line) and the galileon field (Q n , bold dashed line) for 
r c = 20 Gpc, function of redshift. 
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Figure 4: The effective equation of state for the expansion history, defined in terms of the Hubble 
parameter by H/H 2 = — 3(l + w to t)/2, for ACDM (solid), and the galileon model with r c = 10 Gpc 
(dotted), 15 Gpc (dash-dotted) and 20 Gpc (dashed). 
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— 1 at late times. Conversely, w tot approaches — 1 today as r c — * oo. While it may seem surprising 
at first sight that w tot > — 1 at late times, since we argued in Sec. 4 that w w < —1 on the branch of 
interest, one must keep in mind that p n < — the galileon contribution to H, being proportional 
to —(1 + w w )p n , is therefore negative. 

To say a few words about the asymptotic behavior, since p n is negative and decreasing, it 
eventually catches up with p\. At this point the Hubble parameter vanishes, and the universe 
starts to contract. Galileon models thus generically predict a late-time contracting phase. This 
could naturally match onto an ekpyrotic [30, 31, 32, 33, 34], New Ekpyrotic [35, 36, 37, 38, 39] or 
cyclic [40, 41] contracting phase, for instance. See [42] for a review of these models. 

Figure 5 shows the effective BD parameter u^ D , introduced in (4.15), for the whole evolution, 
with r c = 10 Gpc (dotted curve), 15 Gpc (dash-dotted) and 20 Gpc (dashed). At early times the 
uj^ d ~ 10 8 is large, and standard cosmology is recovered. At late times, ~ 1, and the theory 
reduces to scalar-tensor gravity. As expected, the scalar-tensor regime is achieved at earlier times 
for smaller values of r c . 
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Figure 5: The effective BD parameter, u;^, with r c = 10 Gpc (dotted), 15 Gpc (dash-dotted) and 
20 Gpc (dashed). 
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5.2 Comparison with DGP cosmology 



How closely does our AD galileon theory come to reproducing the cosmology of the 5D DGP 
model? To make the comparison it is convenient to think of the H/r c correction term in the DGP 
Friedmann equation (4.8) as an effective energy density component: 

£52*=-— (5 3) 

where we have chosen the normal branch. Figure 6 compares the evolution of |pdgp| i n DGP 
cosmology with the galileon energy density \p n \ in our galileon theory, in each case with r c = 
15 Gpc. For simplicity, we fix the matter, radiation and cosmological constant contributions. We 
see that pdgp and agree remarkably well for the entire expansion history until z — 0, and begin 
to diverge in the future (z < 0) when \tt\ ~ M P ]. 




1+z 



Figure 6: Comparison of the galileon energy density (\p n \) with the effective energy density coming 
from the modification to the Friedmann equation in the DGP model (IpdgpD- We have used 
r c = 15 Gpc in each case and kept the matter, radiation and cosmological constant densities fixed. 
The densities p n and pdgp agree well for the entire expansion history until today, and begin to 
diverge in the future when 7r ~ Mpi. 
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6 Growth of Density Perturbations 



By virtue of its non-minimal coupling to gravity, the galileon enhances the gravitational attraction 
between particles. Translated to the cosmological context, we expect more rapid growth of density 
perturbations. Furthermore, from the discussion of Sec. 2, the galileon enhancement should be 
suppressed at early times, due to self-screening, but should become important once the matter 
density has dropped sufficiently 

Since the matter action is independent of the galileon in Jordan frame, the evolution of matter 
density perturbations, 5 m = 5p m / p m , on sub-Hubble scales is governed by the standard expression 

5 m + 2H5 m = , (6.1) 

cr 

where $ is the Newtonian potential in Jordan frame. The effects of n are all encoded in its 
contribution to the Poisson equation. 

Consider expanding the galileon around its cosmological profile, ir(x, t) = 7t(t) + p(x, t). In the 
Newtonian approximation, we are justified in neglecting time derivatives of tp relative to spatial 
gradients: \ip\ <C |Vy?|. We can thus expand (3.2) to linear order in p as follows, keeping in mind 
that e 2?r / Mpi PS 1 for most of the expansion history, 

-4r 2 (# + 2m) = M^SR . (6.2) 
Meanwhile, from the trace of (3.3), 

-M^6Rfa6T + 6M P] ^-^-. (6.3) 



0/ 



Note that we have neglected the 7r-dependent terms in the second line of (3.3), since the backre- 
action of the galileon is negligible for all times, except in the recent past. Since 5T = — p m S m for 
non-relativistic sources, combining (6.2) and (6.3) gives 



VV 4ttG 2 



M P1 3(3 

where 



aTp 5 m , (6.4) 



2r 2 

Moreover, using (6.4), it is straightforward to obtain the modified evolution equation for density 
perturbations: 

L + 2H5 m = 4irG Pm 5 m (l + ±- V (6.6) 



3[3 
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This Poisson equation for ip exhibits the expected features: at early times, when r 2 \n+2H7i\/Mp\ 3> 
1, the coupling to S m is much suppressed since (5 — > oo. This is the cosmological self-screening 
mechanism, which decouples the galileon from matter inhomogeneities. At late times, on the other 
hand, f3 becomes of order unity, and the galileon couples to matter with strength comparable to 
that of gravity. The evolution equation (6.6) takes on a standard form in the late-time regime, 
except for the fact that the gravitational attraction is enhanced by the galileon factor. 

For comparison, the linearized perturbation equation derived for the full DGP model, derived 
in [43] for spherical top-hat perturbations, takes on an identical form to (6.6), with (3 in this case 
given by 

3DGP l + 2r c H + r 2 r H 2 



(6.7) 



1 + r c H 

(Note that the choice of sign is consistent with the normal branch of DGP.) Clearly the asymptotic 
behavior of (3 DGP for Hr c <C 1 and Hr c 3> 1 agrees with that of (6.5). Figure 7 compares 1/3/? in 
our model with 1/3/9dgp m DGP cosmology, by solving numerically for the respective cosmological 
backgrounds. As in the discussion of Sec. 5.2, we take r c = 15 Gpc in each case and fix the energy 
density in matter, radiation and cosmological term. Given the close agreement shown in the figure, 
we expect our galileon theory to make nearly identical predictions for structure formation as the 
full DGP model. We leave for future work a detailed study of the evolution of perturbations in 
the covariant galileon theory. 
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Figure 7: Comparison of the perturbation enhancement factor 1/3/5 in our galileon model (dotted) 
with the corresponding factor 1/3/3dgp in DGP cosmology (dashed). As in Fig. 6, we have used 
r c = 15 Gpc in each case and kept the matter, radiation and cosmological constant densities fixed. 
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7 Observational Constraints 



A rigorous comparison of the galileon cosmological scenario with observations requires a full likeli- 
hood analysis, which is beyond the scope of this paper. For the purpose of this Section, we restrict 
ourselves to a comparison with a few observables and derive a lower bound on r c . Since r c is the 
only parameter that we allow to vary, our estimate is likely conservative. 



7.1 Mass estimates 

Because of the non-minimal coupling of the galileon, the fractional matter density, fl m , depends 
on 7r, as seen from (5.2): 



P 



e 27r/M Pl 



m ' [ ' 

Therefore, for fixed matter density today, Qm\ the matter density in the past will differ from the 
standard gravity prediction by 

*->" _ e -27T /Mpi ^ 72 ) 



where tiq is the present value of the galileon. And since 7r < on the stable branch, Q m is larger at 
early times than predicted by standard gravity, again keeping fixed. Hence ttq is constrained 
by estimates of the matter density at various redshifts, such as from cluster counts, Lyman-a 
forest and weak lensing observations. (Similar considerations apply to coupled dark matter-dark 
energy models [44].) 

Although the constraints from these observables should be revisited in the presence of a 
galileon, it has been argued that the allowed range of Q$ is generally insensitive to the specifics 
of dark energy [45]. A general analysis combining SNIa Gold data set [46], Wilkinson Anisotropy 
Microwave Probe (WMAP) power spectra [47], and Two-Degree Field (2dF) galaxy survey [48] 
obtained 0.23 < f2m' ) ^ 0.33 [45]. (See also [49, 50].) In the absence of a full likelihood analysis, 
we can obtain a conservative bound on n by requiring that e~ 27r °/ Mpi < 0.33/0.23, or 

|7T | <0.18M Pl . (7.3) 

Figure 8 shows \ir \ for various values of r c , as obtained numerically. We can read off that the 
above bound on | vr 1 is satisfied for r c > 10 Gpc, or roughly 3 times larger than the radius of the 
observable universe. We will see shortly that luminosity distance observations tighten the bound 
to r c > 15 Gpc. 
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Figure 8: Value of the galileon field at the present time as a function of r c . For each value of 
r c , the cosmological term is adjusted to keep Q$ = 0.26 fixed. The dashed line represents the 
threshold value of 0.18 Mpi from mass estimates — see the bound in (7.3). 

7.2 Cosmological Distances 

Next we turn our attention to cosmological distance tests, in particular the luminosity distance 
relation, 



as constrained by Type la supernovae (SNIa). Figure 9 shows the luminosity distance relation 
with = 0.26 for the ACDM model, and for our galileon model with r c = 10, 15 and 20 Gpc. 
In the galileon examples, A was adjusted in order to keep fixed. 

Figure 10 shows the percentage difference between the various galileon examples and the ACDM 
fiducial model. The uncertainties in present SNIa data constrain the luminosity distance to no 
better than ~ 7% over the range < z < 1.5. Therefore, from the percentage differences shown 
in the figure, we see that luminosity-distance observations tighten the bound on r c to 




(7.4) 



r c > 15 Mpc . 



(7.5) 



This constitutes our main constraint on the scale of the modification. 
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Figure 9: Luminosity distance (d L ) as a function of redshift for ACDM (solid), and our galileon 
model with r c = 10 Gpc (dotted), r c = 15 Gpc (dash-dotted) and r c = 20 Gpc (dashed). In each 
case, we have fixed = 0.26. 
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Figure 10: Percentage difference between the galileon examples and ACDM. 



23 



Another distance constraint comes from the angular-diameter distance to the last scattering 
surface, 

1 f Zrcc dz' 

di(z ~ ) = TT^i„ wr (7 - 6) 

which determines the position of the cosmic microwave background (CMB) acoustic peaks. For 
r c > 15 Gpc, the difference in d\(z iec ) compared to ACDM is less than 10%, again keeping 
fixed, which is within current CMB uncertainties. Note that a similar constraint (r c > 3 — S.Sif^ 1 ) 
was obtained recently by confronting the DGP normal-branch cosmology against CMB, SNIa and 
Hubble constant observations [51]. 



8 Conclusions 

In this paper we have studied the cosmology of a galileon field theory, obtained by covariantizing 
the 7r-Lagrangian of the DGP model. Despite being a local theory in 3+1 dimensions, the resulting 
cosmological evolution is remarkably similar to that of the full 4+1-dimensional DGP framework, 
at least for Hr c > 1 (or | -zr j < Mpi). The similarity holds for both the expansion history (Fig. 6) 
and the evolution of density perturbations (Fig. 7). 

In particular, as in the DGP model our covariant galileon theory yields two branches of solu- 
tions, depending on the sign of tt. Perturbations are stable on one branch and ghost-like on the 
other. The effective equation of state for the galileon is phantom-like (w n < —1) on the stable 
branch, and standard (w n > —1) on the unstable branch. A key difference, however, is that the 
galileon field theory does not generate a self-accelerated solution — as shown in Sec. 4.1, the 
would-be self-accelerated solution with H ~ l/r c is spoiled by l/Hr c corrections in our theory. 

An interesting effect uncovered in our analysis is a cosmological version of the self-screening 
(or Vainshtein) mechanism. At early times, Hr c ^> 1, the evolution of n is dominated by the 
self-interactions terms. In turn, this results in n being a negligible component, with its energy 
density suppressed by a factor 0(1/ Hr c ) compared to matter and radiation. This cosmological 
self-screening is crucial in the recovery of standard cosmology at early times. 

Once the expansion rate drops to Hr c ~ 1, however, the galileon becomes an important player 
in the Friedmann equation. A preliminary analysis of observational constraints in Sec. 7 shows that 
the modifications to the expansion history are consistent with the observed late-time cosmology 
provided r c > 15 Gpc. We should emphasize that this bound is most certainly conservative and 
will be relaxed by allowing other cosmological parameters to vary. 

Our analysis offers a host of interesting avenues to explore: 

• Following up on the discussion of the last paragraph, a thorough comparison with observa- 
tions requires a full likelihood analysis, allowing various cosmological parameters (such as 
Vim , h, etc.) to vary. 
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• A study of the implications for structure formation requires a more rigorous treatment of 
inhomogeneities in the presence of the galileon. To probe the non-linear regime, it should 
be straightforward to generalize the N-body simulations of [52, 53] to include the galileon 
field theory studied here. 

• In this work we have focused exclusively on the cubic interaction term that arises in the de- 
coupling limit of DGP. Recently, [11] derived the most general galileon field theory (without 
gravity) which is invariant under the Galilean shift symmetry d^n — > d^n + c M , and whose 
equations of motion are second order. It was shown in [54] that the equations remain of 
second order in the presence of gravity, provided n is suitably coupled to gravity. It would 
be very interesting to extend our analysis to include these higher-order interactions terms. 

• The extension of our 4D theory to include multiple galileon fields should offer a reliable proxy 
for the cosmology of higher-dimensional DGP models, such as Cascading Gravity [55, 56, 57]. 
In this construction, our 3-brane lies within a succession of higher-dimensional DGP branes, 
embedded in one another within a flat bulk space-time. The corresponding AD covariant 
theory should therefore include 2 interacting galileon fields, each with its own r c scale. 
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